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ABSTRACT 



To understand the formation and evolution of solar-type stars in the solar neighborhood, we need to measure their stellar parameters to high 
accuracy. We present a catalogue of accurate stellar parameters for 45 1 stars that represent the HARPS Guaranteed Time Observations (GTO) 
"high precision" sample. Spectroscopic stellar parameters were measured using high signal-to-noise (S/N) spectra acquired with the HARPS 
spectrograph. The spectroscopic analysis was completed assuming UTE with a grid of Kurucz atmosphere models and the recent ARES code 
for measuring line equivalent widths. We show that our results agree well with those ones presented in the literature (for stars in common). 
We present a useful calibration for the effective temperature as a function of the index color B-V and [Fe/H]. We use our results to study the 
metallicity-planet correlation, namely for very low mass planets. The results presented here suggest that in contrast to their jovian couterparts, 
neptune-like planets do not form preferentially around metal-rich stars. The ratio of jupiter-to-neptunes is also an increasing function of stellar 
metallicity. These results are discussed in the context of the core-accretion model for planet formation. 



Key words. Methods: data analysis - Technics: Spectroscopy - stars: abundances - 
stars: planetary systems: formation - galaxy: chemical evolution 
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1. Introduction 

Extra-solar planets have continued to be discovered, since 
the first detection of a gi ant planet orbiting a solar-like 
star ( Mayor & Oueloz 1995|). including several planets in the 
Neptune-ma ss regime (e .g. McArthuret al. 2004; 'Santos et al.' 

12004; Vogt et al, 2005- Rivera et al. 2005- 



date, which has achieved an instrumental radial-velocity ac- 
curacy superior t o 1 ms"' for many years (.Mayor et al.i,2003l: 
Lovis et al.ll2005l) . 



2004b; Buder et al 



Bonfils et al 



2005aL 120071 lUdrv et alJl2006f 120071: iLovis et al 



2006t iMelo et alj l2007 l). This has been possible because of 



to the development of a new generation of instruments capa- 
ble of radial-velocity measurements of unprecedented quality. 
One such instrument is undoubtedly the ESO high-resolution 
HARPS fiber-fed echelle spectrograph, especially designed 
for planet-search programmes and asteroseismology. HARPS 
has proved to be the most precise spectro-velocimeter to 



Send offprint requests to: S. G. Sousa: sousasag@astro.up.pt 

* Based on observations collected at La Silla Observatory, ESO, 
Chile, with the HARPS spectrograph at the 3.6-m telescope (072. C- 
0488(E)) 



A by-product of these planet searches are the high-quality 
spectra, of high S/N and high-resolution, that are collected in 
every run in a consistent and systematic way for selected can- 
didate stars. The result is a vast quantity of spectra, mainly of 
solar-type stars, which can be analyzed in a consistent way for 
the determination of accurate spectroscopic stellar parameters. 

We present a catalogue of spectroscopic stellar parameters 
for the HARPS Guaranteed Time Observations (GTO) "high 
precision" sample. In Sect. 2 we describe the observations with 
HARPS. Section 3 describes the procedure used to derive ac- 
curate spectroscopic stellar parameters. In Sect. 4, we present 
the results for our sample of 45 1 solar-type stars. We then com- 
pare, in Sect. 5, our results with those of several other authors 
obtained using different methods. In Sect. 6, we present a use- 
ful calibration of the effective temperature as a function of B-V 
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and [Fe/H]. Finally, we perform a simple statistical analysis of 
our results in the context of the metallicity-planet correlation. 



2. Stellar Sample & Observations 

The objective of the HARPS "high-precision" GTO project is 
to detect very low-mass extra-solar planets by increasing the 
radial velocity measurement accuracy below 1 ' . The stars 
in this project were selected from a volume-limited stellar- 
sample obse rved by the COR ALIE spectrograph at La Silla 
observatory ( Udry et aDl2000f) . These stars are slowly-rotating, 
non-evolved, and low-activity stars that presented no obvious 
radial-velocity variations at the level of the CORALIE mea- 
surement precision (typically 5-10 mi"' and dominated by 
photon-noise). Southern planet hosts were also added to this 
program to search for further low-mass planets. Therefore, a 
significant number of stars (66) from this catalogue are already 
confirmed planet hosts. For mor e details we point t he reader 
to a description of this sample bv lMavor et alj (I2003I) . In sum- 
mary, the sample is composed of 45 1 FGK stars with apparent 
magnitudes that range from 3.5 to 10.2 and have distances less 
than 56 parsec. 

Asteroseismology work using HARPS has proved that the 
precision of the instrument is no longer set by instrumental 
characteristics but rather b y the stars themselves ( Mayor et al.l 
)03t iBouchv et al.ll2005l) . Even a "quiet" G dwarf shows os- 
cillation modes of several tens of cmi ' each, which might 
add up to radial-velo city amplitudes as large as several mi"' 



Table 1. Examples of atomic parameters and measured solar 
equivalent widths for Fe n and Fe i lines. 



(IBouchv et alJ l2005l) . As a consequence, any exposure with 
a shorter integration time than the oscillation period of the 
star, might fall arbitrarily on a peak or on a valley of 
these mode interferences and thus introduce additional radial- 
velocity "noise". This phenomenon could seriously compro- 
mise the ability to detect very low-mass planets around solar- 
type stars by means of the radial-velocity technique. 

To average out stellar oscillations, the observations are de- 
signed to last at least 15 minutes on the target, splitting them 
into several exposures as required, to avoid CCD saturation. 
The final result is a high-quality spectrum for each star in the 
sample. 

The individual spectra of each star were reduced using the 
HARPS pipeline and combined using IRAF after correcting 
for its radial velocity. The final spectra have a resolution of R 
~ 110000 and signal-to-noise ratios that vary from ~ 70 to 
~ 2000, depending on the amount and quality of the original 
spectra (90% of the spectra have S/N higher than 200). 



3. Stellar parameters - Spectroscopic Analysis 

We followed the procedure presented in I Santos et alJ (l2004bh 
for deriving stellar parameters and metallicities, where [Fe/H] 
is used as a proxy for overall stellar metallicity. As reported in 



' IRAF is distributed by National Optical Astronomy 
Observatories, operated by the Association of Universities for 
Research in Astronomy, Inc., under contract with the National 
Science Foundation, U.S.A. 
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Sousa et alj (120071) however, we attempted to increase the num- 
ber of spectral Unes and decrease the errors involved in the de- 
termination of stellar paramete rs. Therefore, we u sed the large, 
extended line-list presented in lSousa et al. I (l2007h and applied 
this method to the sample. 

The spectroscopic analysis was completed assumi ng LTE, 
and using the 2002 version of the code MOOC0 ( SnedenI 
I973h and a grid of Kur ucz Atlas 9 plane-parallel model at- 
mospheres ( Kurucziri993 ). Stellar p arameters and metal l icities 
were derived as in previous works dSantos et al. I l2004b[l2005l 



Sousa et al. I l2006h . based on the equivalent widths of Fe i and 
Fe II weak lines, by imposing excitation and ionization equiUb- 
rium. 

3.1. A large list of stable lines 

We cons idered the Une list used for the spectroscopic analysis 
given in lSousa et alj (l2007h . As a first step, to check the qual- 
ity of the line list, we derived preliminary stellar parameters 
for our sample. The stability of the lines were then checked 
by comparing the abundances of each individual Une with the 
mean abundance obtained for each star. 

Since the stars in our sample have a significant range of 
spectral type with efifective temperatures ranging from -4500 
K to ~6400 K, we expect a significant different local envi- 
ronment around each line, due to variations in blending ef- 
fects. To take this problem into account, we divided our sample 
into three subsamples, taking into consideration the prelimi- 
nary determination of the effective temperature (Teff < 5300^", 
5300/: < Teff < 6000/:, Teff > 6000/:). 

This allows the search for systematic non-stable lines in the 
large list for the three ranges of stellar spectral class. The lines 
that showed a systematic offset from the mean value, that was 
larger than 1.5 sigma were considered "bad" lines, which pro- 
vide false abundance measurements either due to line-blending 
events or because the continuum could not be fitted readily in 
a particular spectral region. Lines at wavelengths lower than 
4500A were removed because, for cooler stars, this region 



^ The source code of MOOG2002 can be downloaded at 
|http://verdi.as.utexas.edu/moog.html] 
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Fig. 1. Sele c tion o f the stable lines from the list provided by 
Sousa et alj (120071) . We present the mean difference between 
each Une abundance and the mean abundance of the respective 
star as a function of the wavelength. We divided the analysis 
into three different ranges of effective temperatures. The dashed 
lines represent the 1 .5 sigma threshold. Lines in the shaded area 
(< 4500A) were also removed since this region was strongly 
affected by blending effects, in particular for cooler stars. 



is highly line crowded and therefore more likely to produce 
blending effects. These lines were removed from the original 
large line list. 

Table [T] shows the final list composed of 263 Fe i and 36 
Fe II weak lines to be used on HARPS spectra for the determi- 
nation of the stellar parameters and iron abundance. 

We note that this procedure is difficult to apply to cooler 
stars, because the spectra of these stars are crowded with lines 
or the atomic parameters that we adopt become inaccurate 
as we study spectral classes increasingly different from solar 
These two effects can be seen clearly in the dispersion of mea- 
surements shown in Fig. [T] (upper panel) that is significantly 
higher than for the two other cases (middle and lower panel). 

Due to these problems, the solutions for cooler stars are 
more difficult to converge into the "true" stellar parameters. In 
a few cases, we had to discard some lines from the final stellar- 



3.2. Using ARES 

The equivalent widths of the lines were measured automatically 
in all spectra using the AREdl code (A utomatic Routine f or 
line Equivalent widths in stellar Spectra - Sousa et al ]l2007h . 

The ARES code considers a one-dimensional spectrum, a 
list of absorption lines to be measured, and some calibration pa- 
rameters. First, ARES estimates the continuum about each in- 
dividual line and calculates a local normalization. The location 
and number of lines for which a Gaussian fit can be attempted, 
are then determined. After fitting the local, normalized spec- 
trum, ARES uses the fit parameters to compute the equivalent 
width. The result is then outputted to a file. 

A number of ARES parameters were fixed during the anal- 
ysis of all spectra in our sample: smoothder - 4, space - 3, 
lineresol — 0.07, miniline — 2. The smoothder value was used 
to reduce the noise in the computed derivatives of the local 
spectrum to find the number and location of the lines to be 
fitted. The space parameter was set to be the total length of 
spectrum about each line required to fit the local continuum; 
therefore, for a value of 3 A, ARES fitted the continuum over a 
total interval of 6 A. The lineresol parameter defines the mini- 
mum separation allowed for consecutive lines. This parameter 
is used to limit the effect of noise on the determination of the 
number and location of lines, i.e. if ARES identify lines that 
are close together according to the value indicated in lineresol, 
then it ignores one line and then measures a mean position. 
The miniline parameter is the lower value of EW that will be 
returned by ARES in the output file. 



Table 2. Dependence of the rejt parameter on the S/N of 
HARPS spectra. ARES users should use this table as a refer- 
ence for the rejt parameter 



S/N condition 


rejt 


S/N condition 


rejt 


S/N < 100 


0.985 






100 < S/N < 125 


0.990 


250 < S/N < 300 


0.995 


125 < S/N < 150 


0.991 


300 < S/N < 500 


0.996 


150 < S/N < 200 


0.992 


500 < S/N < 800 


0.997 


200 < S/N < 225 


0.993 


800 < S/N < 1500 


0.998 


225 < S/N < 250 


0.994 


1500 < S/N 


0.999 



As described in ISousa et al.l (l2007h . the most important 
parameter for the correct automatic determination of equiva- 
lent widths using ARES is the rejt parameter. This parameter 
strongly depends on the signal-to-noise ratio (S/N) of the spec- 
tra. Values close to 1 are indicative of lower spectrum noise 
levels. Table |2] shows how the parameter was defined by con- 
sidering the spectrum S/N (measured at 6050 A). Future ARES 



^ The ARES code can be 
pittp://www.astro.up.pt/~sousasag/ares| 



downloaded 



at 
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Table 3. Sample table of the HARPS GTO "high precision" spectroscopic catalogue. We provide the star name, effective tem- 
perature, surface gravity, micro turbulence, [Fe/H], the number of iron lines used in the spectroscopic analysis, a mass estimate, 
the surface gravity computed from the parallaxes, the luminosity of the star, and an indication of whether it hosts a planet. In the 
electronic table, we also indicate the error in the luminosity and the computed absolute magnitude of the stars. 



Star ID 


Tefif 




log 


8 spec 






[Fe/H] 


N(Fe i,Fe ii) 


Mass 


^Og ghipp 


Luminosity 


planet host? 




[K] 




[cms ^] 


[kms-i] 






[Mq] 


[cm s"^] 


ILq] 




HD 63454 


4840 ± 


66 


4.30 


+ 0.16 


0.81 


± 0.14 


0.06 ± 0.03 


227, 27 


0.66 


4.55 


0.25 


yes 


HD 63765 


5432 ± 


19 


4.42 


+ 0.03 


0.82 


±0.03 


-0.16 ±0.01 


249, 32 


0.83 


4.51 


0.55 


no 


HD 65216 


5612 ± 


16 


4.44 


±0.02 


0.78 


± 0.03 


-0.17 ±0.01 


256, 34 


0.87 


4.48 


0.71 


yes 


HD 65277 


4802 ± 


88 


4.43 


±0.18 


0.55 


±0.34 


-0.31 ±0.04 


257, 32 


0.49 


4.49 


0.21 


no 


HD 65562 


5076 ± 


47 


4.39 


±0.09 


0.45 


±0.18 


-0.10 ±0.03 


252, 34 


0.72 


4.51 


0.37 


no 


HD 65907A 


5945 ± 


16 


4.52 


±0.02 


1.05 


±0.02 


-0.31 ±0.01 


256, 34 


0.93 


4.37 


1.24 


no 


HD 66221 


5635 ± 


25 


4.40 


±0.04 


0.92 


±0.03 


0.17 ±0.02 


260, 35 


0.96 


4.40 


0.95 


no 


HD 66428 


5705 ± 


27 


4.31 


±0.06 


0.96 


±0.03 


0.25 ± 0.02 


258, 35 


1.01 


4.31 


1.28 


yes 


HD 67458 


5891 ± 


12 


4.53 


±0.02 


1.04 


±0.01 


-0.16 ±0.01 


256, 35 


0.98 


4.45 


1.02 


no 


HD 68607 


5215 ± 


45 


4.41 


±0.08 


0.82 


±0.08 


0.07 ± 0.03 


251,35 


0.81 


4.52 


0.45 


no 


HD 68978A 


5965 ± 


22 


4.48 


±0.02 


1.09 


± 0.02 


0.05 ± 0.02 


260, 34 


1.08 


4.43 


1.24 


no 


HD 69655 


5961 ± 


12 


4.44 


±0.03 


1.15 


±0.01 


-0.19 ±0.01 


249, 34 


0.98 


4.36 


1.33 


no 


HD 69830 


5402 ± 


28 


4.40 


±0.04 


0.80 


±0.04 


-0.06 ± 0.02 


255, 36 


0.82 


4.46 


0.60 


yes 


HD 70642 


5668 ± 


22 


4.40 


±0.04 


0.82 


±0.03 


0.18 ± 0.02 


253, 36 


0.98 


4.42 


0.95 


yes 



users may use this table as a reference for the most appropri- 
ate choice of the rejt parameter. This table was compiled for 
HARPS spectra, but the data should be applicable to other high- 
resolution spectra. 

3.3. Stellar mass and luminosity 



Stell ar masses were 
e.g. ISantos et"al 



estimated, as in previous works 

(l2004bl). by i nterpo lating the theoreti- 



cal isochrones of ISchaller et al.l (Il992h and ISchaerer et al 



(11993 b a), using My computed for the Hipparcos parallaxes 
and V magnit udes (lESAlll997h . a bolometric correction from 
lFloweil(ll996h . and the effective temperature derived from the 
spectroscopic analysis. We adopt a typical relative error of 0. 10 
Mq for the masses. In some cases, the value of mass was not 
determined because the calculation would have required a large 
extrapolation of the existing isochrones. 

The luminosity was computed by considering the 
Hipparcos parallaxes, V magnitude and the bolometric coiTec- 
tion. Its error was derived based on the parallax eiTors from 
Hipparcos, which was the main source of uncertainty in the cal- 
culation of luminosity. The eiTor in the luminosity is available 
in the electronic table only. 



4. The spectroscopic catalogue of the HARPS GTO 
"high precision" sample 

Figure |2] presents some characteristics of the sample. The top 
plot shows the distribution of the sample on the Hertzsprung- 
Russell diagram, where we represent evolutionary tra cks for 
1.0, 1 1, and 1.2 comp uted using the CESAM code (iMorell 
1997t Marques et al.ll2007h K Moreover, we present the typical 



error boxes on this specific diagram, where the error in the lu- 
minosity is dependet mostly on the error in the parallax of each 
star, and the error in the temperature is derived from our spec- 
troscopic method. This plot shows that the sample is composed 
mainly of main-sequence solar-type stars. In the bottom plot, 
we present the metallicity distribution that has a mean value of 
about -0.09 dex. This is co mpatible with t he com parison sam- 
ples presented in the work of lSantos et al. I (l2004bh . when taken 
into consideration the amount of planet hosts (more metal-rich 
on average) that were added to the original HARPS GTO "high 
precision" sample. 

Table |3] shows a sample of the catalogue with some of the 
stellar parameters that have been determined. This catalog is 
electronically available on CDS and also at the author web- 
pag^ 



5. Comparison with previous works 

We compared our results with those of other authors to 
acess the consistency between different methods, using spec- 
troscopy or photometry to deriv e stellar paramet e rs. W e 
used the spe c trosco p ic results of Edyardsson et alj ( 1993b. 
Bensbvetal. (20031). IValenti & Fischeil (l2005b . ISantos et al. 
(2004b). Fuhrmanni (20041) . and the photometric results of 



Nordstrom et al. ( 20041) . We also checked the consistency with 



results derived using the Infr a-Red Flux M ethod (IRFM), fol- 
lowing the work of Casagrande et al.l (120061) . 



From this point on, all differences presented are relative to 
our own work, i.e. "other works"-"this work". 



http://www.astro.up.pt/corot/models/ 



^ http://www.astro.up.pt/~sousasag/harps_gto_catalogue.html 
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Fig. 3. Comparison of our spectroscopic results for the effective temperature with other values found in the literature. We also 
show the comparison with IRFM using either the Kurucz or Phoenix models. 



5.1. Effective Temperature 

For the effective temperature the results of the comparison, pre- 
sented in Fig. ID are compatible with other spectroscopic stud- 
ies. We found a mean difference o f only -18+ 47 K for the 20 



stars in common with the work of iBensby et al.l (l2003h . These 
authors used the same procedure as this paper for determining 
the effective temperature, but in their case for different stellar- 
atmosphere models (Uppsala MARCS code). In their work, 
they measured equivalent widths from spectra taken with the 
FEROS spectrograph on the 1.52 m ESO telescope. 

[Fuhr mann (2004) used a different procedure for the deter- 
mination of the effective temperature, based on hydrogen line 
profiles. However, our and their results are in excellent agree- 
ment for the 15 stars that we have in common, for which we 
fin d a mean difference of -h1 ± 44 K. 

Valenti & Fischer (l2005h presented stellar parameters for 



1040 F, G, and K dwarfs. We were able to complete a com- 
parison using data for a larger number of stars and found good 
agreement with a mean difference -16 + 44 K in the effective 
temperature. However, a small offset can be seen in the plot for 
higher temperatures (Teff > 6000 K). These authors used spec- 



tral synthesis for data taken with HIRES, UCLES, and the Lick 
Observatory Hamilton Echelle Spectrometer. 

We compared our measurements derived here with our own 
previous results ( Santos e t al. 2004b), which were obtained us- 
ing the same procedure but a different line list. The compari- 
son indicated good agreement for effective temperatures with a 
mean difference of -18 +57 K. For this specific comparison , 
most stars in common with the sample of lSantos et alj (l2004bl) 
were observed with the CORALIE spectrograph and only a few 
with the FEROS spectrograph. 

In their careful spectroscopic analysis of 189 nearby field F 



and G disk dwarfs, Edvardsson et al.l (Il993h determined effec- 
tive temperatures using photometry. For the stars in common, 
we found very similar effective temperatures, with a mean dif- 
ference of only -27 + 71 K. 

We also compared our effective temperatures with those 
in Nordstrom et al. (2004) that used the IRFM-based calibra- 
tion of lAlonso et al.l (Il996l) . T here is a systematic devi ation be- 
tween our results and those of iNordstrom et al.l ( 120041) . with a 
mean difference of -108 + 68 K. In particular, a clear offset 
appears to be present for effective temperatures above 5500 K. 
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Fig. 4. Comparison of our spectroscopic results for surface gravity with other measurements in the literature. 



We do not know th e reason for this, but it is interesting to note 
that the lAlonso et al. ( 1996) temperat ure scale is based on the 
Vega absolute calibration derived by lAlonso et alj ( Il994l) . In 
that work, it was found that the Vega absolute calibration de- 
rived from hot (^eff > 6000 K) and cool (T^s < 5000 K) stars 
differed by a few percent, which implied that effective temper- 
atures for cool and hot stars were not on exactly the same scale. 
The Vega absolute calibration used by Alonso et al. (1996) was 
a weighted average of that for both hot and cool stars: it is in 
fact interesting to note that the offset in Fig. [3] appears approx- 
imately where the method showed some problems. 



5.2. Effective Temperature - The IRFM 

Casagrande et al. (2006) implemented a revised form of the 
IRFM, as described in their paper. They found good agreement 
with various spectroscopic temperature scales; they explained 
a ~ 100 K systematic difference with respect to previous im- 
plementation of the IRFM, as being most likely due the adop- 
tion of t he new Vega zero points and absolute calibration of 
2MASS dCohen et alj|2003) . The comparison with a large set 
of solar analogs confirm that the temperature sca le is calibrated 
well for solar-type stars (ICasagrande et al ] l2006l) . 

The main concept behind the IRFM is to assume that the ra- 
tio between the bolometric and infrared monochromatic fluxes 
is a sensitive indicator of effective temperature. The bolomet- 
ric flux is recovered using multi-band photometry, right across 
the optical and near-infrared. The flux outside the photometric 
filters is estimated using model atmospheres, and the infrared 
monochromatic fluxes are computed from 2MASS JHR^ pho- 
tometry. Accurate infrared photometry is crucial for recover- 
ing effective temperatures precisely and only stars with total 



photometric eiTors (i.e. "j_"H-"h_"H-"k_msigcom" as given from 
2MASS) smaller than 0. 10 mag are used. 

For 66 stars in the HARPS GTO sample we obtained 
accurate Johnson-Cousins BV(RI)c colors from t he Ge neral 
Catalogue of Photometric Data (Mermilliod et aL| |l997|) and 
JHKs from 2 MASS, and executed the IRFM as described in 



Casagrande et al 



(l2006h . To increase the number of stars for 
which it was possible to apply the IRFM, we tested the results 
only if Johnson BV and 2MASS JHKs magnitudes were used. 
We found the average difference to be AT,,// = 5 + 15 K, im- 
plying that the IRFM could be applied to aU H ARPS GTO stars 
which with Hipparcos Johnson photometry (lESAi(l997l) . We 
also checked that if Hipparcos Tycho BtVt magnitudes, trans- 
formed onto the Johnson system (Bessell 2000), were used in- 
stead that the differences in our results would be negligible. 

Figure [3] (three right plots) shows the comparison with ef- 
fective temperatures derived by the spectroscopic analysis de- 
scribed in Sect. 3. The agreement is excellent, corresponding to 
a mean difference of - 1 + 63 K. We measure a slight disagree- 
ment when we consider cooler stars alone, which may due to 
inconsistencies in stellar model atmospheres at the lowest tem- 
peratures covered by the present study. In the case of t he IRFM, 
we use the ATLAS9-ODFNEW models (Castelli & Kurucz 
I2OO3I) to recover the bolometric flux that is missing from our 
multi-band phot ometry. We also tested th e use of the latest 
Phoenix models dBrott & Hauschildl l2005h . which include an 
extensive molecular-line list, which is particularly important 
for cooler stars. Using the IRFM, we found that results derived 
using the ATLAS9-ODFNEW models agree incredibly well 
with those for the Phoenix models, corresponding to a mean 
difference in temperature of only -5 ± 11 K. Such agreement 
is in part due to our implementation of the IRFM, which uses 
more observational data and highlights the high degree of con- 
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Fig. 2. In the top panel, we show the distribution of the sam- 
ple stars in the H-R diagram. The filled circles represent the 
planet hosts in our sample. We also plot some evolutionary 
tracks computed with CESAM for a 1.0, 1.1 and 1.2 Mq. In 
the bottom panel, we present the metallicity distribution of the 
sample. 



sistency between physical descriptions in independent spectral 
synthesis codes. 

The comparison with the IRFM in Fig. [3] implies that our 
spectroscopic method decreases its accuracy when applied to 
stars cooler than ~ 4800 K. This conclusion is expected since 
the method described in Sect. 3 is based on a differential anal- 
ysis with respect to the Sun; it is also supported by considering 
the error determination for the stars: since errors are based on 
the distribution of the abundance determination, the dispersion 
is larger for cooler stars. Above 4800 K, we find, however, very 
good agreement. 

5.3. Surface Gravity 

Figure |4] shows the comparison between surface gravity mea- 
surements derived in our work and those from other au- 
thors. We also computed spectroscopy surface gravities using 
Hipparcos para llaxes and th e estimated masses of stars as it 
is described in ISantos et al . (2004b). For the comparison be- 
tween these computed gravities and the spectroscopic gravities 
it is observed a dispersion and offset for lower gravities. This 
is probably due to uncertainties in the estimated value of stellar 
mass, which is required in this approach to compute the gravity. 
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16 
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Fig. 5. Comparison of our spectroscopic measurements of 
[Fe/H] to others found in the literature. 



The large mass uncertainty can be attributed to the considerable 
number of degeneracies in stellar interior models: it is possible, 
for example, to describe the data for a star in the H-R diagram 
using several models with di fferent initial conditions (see e.g. 
Fernandes & Monteiroll2003h . It is interesting that some high 
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gravity values (close to 5.0 dex) are computed using the paral- 
laxes. These high values are not expected, since we are dealing 
with dwarf main sequence stars. 



The surface gravities presented by iBensbv et al.l (l2003h 



were determined using the stellar mass and par allax of stars 
In Fig , m we compare the values of log g from iBensby et al 
(120031) with our surface gravities computed from parallaxes. 
We observe perfect agreement with a mean difference of 
-0.02 + 0.05 dex. 

We also find an excellent agreement with the surface grav- 

I IT ' 

ities determined by Fuhrmann (2004) in addition t o our old re- 



sults u sing the same method but a smaller line list (iSantos et al 
2004bl). 



lEdvardss on et alj (119931) measured the surface gravity us- 
ing a different method, which used the Balmer discontinuity 
index. When comparing their values with our spectroscopic 
surface gravities, we find a rather large deviation with a mean 
difference of -0. 14 + 0. 11 dex. 

The mean d if ferenc e between values obtained by 
Valenti & Fischen (l2005h . using spectral synthesis, and 
ourselves is only 0.04 + 0.13 dex, which indicates good 
agreement between the measurements. There is, however, 
significant dispersion in this comparison. For a few stars. 



Valenti & Fischen (l2005h derive unusually high values of grav- 
ity; in the lower regime of gravity, they also measure higher 
values than the o nes we present i n this work. A comparison 
between Va lenti & Fischei (l2005h gravities and independent 
computation of gravity using Hipparcos parallaxes reveals a 
larger offset for smaller gravities. We note that this does not 
imply that gravities given by these authors are incorrect, but 
that the offset is greater than for our values when comparion to 
the computed gravities using parallaxes. We also recall that the 
mass estimates have a large uncertainty, which may contribute 
to the observed offsets. 



5.4. Metallicity 

Figure |5] shows the comparison of the metallicity estimates in 
this work with those derived in other works. As for the other 
parameters, all metallicities agree w ell within the errors. The 
results from tedvardsson et al. ( 1993 ) differ the most from our 
values (with only 0.05 dex), in particular for metallicities above 
-0.1 dex. 



6. A calibration of the effective temperature as a 
function of B-V and [Fe/H] 

Using the parameters presented in Table [3] and the B-V value 
taken from the Hipparcos catalogue (ESA| |1997|) . we derived 
a new calibration of the effective temperature as a function of 
B-V and [Fe/H]. The result is illustrated in Fig. |6] and the cali- 
bration is expressed by; 



Teff = 9114 - 6827(B-y) + 263?.(B-Vy + 368[Fe/H]. (1) 

The standard deviation of the fit is only 47 K, illustrating 
the quality of the relation. Such a result was already observed in 



4000 
0, 




[Fe/H] = 0.5 
[Fe/Hl = 0.0 
[Fe/H] = -0.5 



0.8 0.9 

(B-V) 



Fig. 6. Calibration of the effective temperature as a function of 
the color index B-V and [Fe/H]. The 3 fitted lines correspond 
to lines with constant values of [Fe/H] (-0.5, 0.0, 0.5). 



a caUbration completed bv lSantos et al. I (l2004bl) . using a much 
smaller number of stars. This calibration can be useful and ap- 
plied to stars without the need for a detailed spectroscopic anal- 
ysis, with the guarantee that the result will be in the same effec- 
tive temperature scale. This calibration is valid in the following 
intervals: 4500/: < T.jy < 6400/:, -0.85 < [Fe/H] < 0.40, 
and 0.51 < B-V < 1.20. The small dispersion also attest to the 
quality of the metallicity values derived in this work. 

7. Planet-hosts in the sample 



The well-established correlation between the prese nce of a 
giant p lanet and the meta llicity of its host star (San tos et al 
2004bi: Ipischer & Valentii l2005b is illustrated by our results 



This trend is observed in Fig. |7] where the metallicity distribu- 
tion of jovian-like planets is clearly located towards the metal- 
rich regime. The subsample of 66 planet hosts has a mean 
metallicity of +0.09 dex. For comparison, the sample presented 
in this work has a mean value of -0.12 dex, if we only include 
stars without planets, and -0.09, if we include all the stars (see 
also Fig.|2]l. 

From the 66 planet hosts belonging to our sample, three 
host neptunian planets (HD4308, HD69830 and HD160691) 



The first two of these only harbor neptunian planets (lUdrv et al 



I2OO6; Lovis et al. 2006), while the third also has Jupiter-like 
planets (iPepe et al.l 120071) . Others Neptune-like planets were 
discovered with HARPS, but do not belong to this sample. In 
FiglT] we present the metallicity distribution of these two types 
of host stars. Al thought the numbers are small, we find a wider 
spread of metallicities for stars hosting Neptunian-like planets. 

In Table |4] we compile a list of [Fe/H] values for all 
known planet hosts orbited by neptunes or super-earths. This 
list was compiled using the extrasolar-planets encyclopaedia 
(http://www.exoplanet.eu), for which we select the discovered 
planets with masses (Mp sin /) less than 25 M^. In the table, we 
indicate whether each star also hosts jovian planet(s) ("yes"), 
or if it only hosts neptunian planets or super-earths ("no"). The 
[Fe/H] values are listed in the table for each star. 
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Fig. 7. Metallicity distribution of jovian planet hosts and nep- 
tunian planet hosts within the HARPS GTO "high precision" 
spectroscopic catalogue. 



To increase the number of neptunian-planet hosts, we must 
include M stars presented in the table. We note that the metal- 
licity determination of M stars can be difficult; we therefore use 
the mean of the [Fe/H] values quoted in the literature. 

We should mention that we do not expect the inclusion 
of M-dwarfs in our comparison to strongly compromise or 
bias our conclusions. First, the metallicity distribution of M- 
dwarfs in the solar neighborhood differs only slightly from that 
of field FGK-dwarfs. The offset on the metal distribution is 
small (~0.07 dex), a nd possibly within the margins of errors 
dBonfils et al .12005 a'). Secondly, the fact that these metallicities 
are more uncertain will produce a scatter, but likely not a sys- 
tematic offset, in the [Fe/H] distribution of stars with Neptune- 
like planets. Finally, it is reasonable to assume that the influ- 
ence of metallicity on planet-formation efficiency is indepen- 
dent of either host mass or spectral type. 

Using the values of [Fe/H] listed in Table |4] we compared 
the metallicity distribution of the spectroscopic catalogue pre- 
sented in this paper, which corresponds to the [Fe/H] distri- 
bution of the solar neighborhood, with two different distribu- 
tions. The first, composed of all stars presented in Table|4] 
and the second, composed only of stars that do not host jo- 
vian planets ("no" in the table). The result can be observed in 
Fig-El The metallicity distribution of the catalogue presented 
in this work has a mean metallicity of < [Fe/H] >- -0.09, 
while the metallicity distribution of stars known to harbor at 
least one neptunian planet corresponds to a mean metallicity of 
< [Fe/H] >- -0.03. Finally, stars known to host only neptu- 
nian planets have a mean metallicity of < [Fe/H] >= -0.21. 

We performed a Kolmogorov-Smirnov test to check the 
probability that samples are part of the same population. 
We compare the sample composed of all stars in TablelH to 
that composed only of giant-planet hosts in the catalogue. 
According to a Kolmogorov-Smirnov test, there is a probability 
of 3% that the two samples belong to the same population. We 
applied the same test to stars that host only neptunian planets 
("no" in the table) and found a probability of only 0.5%. This 



Table 4. Neptunian hosts and their metal content. We also list 
the mass (in earth masses) of the least massive planet orbiting 
the star, as well as an indication of whether there is one or more 
jovian companions. 



star 


[Fe/H] 


jov.? 


Mp sin ( 


[Fe/H] reference 


G1581 


-0.33 


no 


5.09 


Bean et al. ("2006") 




-0.25 






Bonfils et al. f2005b) 




-0.29 






Mean value 


Gl 876 


-0.12 


yes 


5.72 


Bean et al. (2006^ 




-0.03 






Bonfils et al. ('2005b) 




-0.08 






Mean value 


HD69830 


-0.06 


no 


10.49 


This work 


HD 160691 


0.30 


yes 


13.99 


This work 


Gl 674 


-0.28 


no 


11.76 


Bonfils et al. ("2007) 


55 Cnc e 


0.33 


yes 


10.81 


Santos et al. ('2004b) 


HD4308 


-0.34 


no 


14.94 


This work 


HD 190360 


0.24 


yes 


18.12 


Santos et al. C2004b) 


HD219828 


0.19 


yes 


20.98 


Melo et al. f2007) 


GJ436 


-0.32 


no 


22.89 


Bean et al. (2006) 




0.02 






Bonfils et al. f2005b) 




-0.15 






Mean value 


Gl 176 


-0.10 


no 


24.16 


Endl et al. (2007) 



result, although preliminary, implies that these samples are un- 
likely to belong to the same population. 

Although the sizes of the datasets are small, we attempt to 
estimate a lower limit value of the frequency of Neptune-like 
planets (or super-earths) found, as a function of stellar metal- 
licity. For each of the three metallicity bins [Fe/H] < -0.15, 
-0.15 < [Fe/H] < 0.15, and [Fe/H] > 0.15, we computed flie 
ratio of the number of stars with Neptune-like planets in Table 
in to the number of stars in our entire sample. In this estimate, 
we assume that the metallicity distribution of the 45 1 stars rep- 
resents the metallicity distribution of the planet search samples 
used to find the planets listed in Table |4] 

The results show that the frequency of stars hosting nep- 
tunians is 1.8%, 2.0%, and 5.1%, respectively for each of the 
metallicity bins mentioned above. These numbers may, how- 
ever, be biased towards the high-metallicity regime, since sev- 
eral stars in Table |4] also have jovian planets. We further note 
that we attempted to find Neptune-like planets orbiting some 
of the FGK stars in Tab leBl because of their high metallic- 
ity (e.g. Melo et aPbOOTh . In other cases, the approach taken 
to detect a Neptune-like planet assumed that th e star hosted 
already at least one jupiter-l ike companion (e.g. ISantos et al 



ioy 

3[ 



200M McArthuret 



ill 



20041) . If we consider stars hosting only 
Neptune-like planets, these frequencies change to 1.8, 1.5, and 
0.0% (there are no stars that host Neptune-like planets, that do 
not host Jupiter-like planets with [Fe/H]> -0.06). 

These numbers imply that the few low-mass planets found 
appear to follow a trend in metallicity that differ from well- 
established rela tions for giant planets first pointed out by 
IUdrvetal](l2006i) . this result, is supported by planet-formation 
model s based on the core-accretion paradigm llda & Linll2004l: 



Benz et al. 2006) . These works suggest that neptunian planets 



should be found in a wider range of stellar metallicities. Lower- 
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Fig. 8. Metallicity distribution for the sample presented in this 
work (dotted line), for stars hosting neptunian planets (dashed 
line), and for stars exclusively hosting neptunian planets (full 
line). The mean metallicity of each distribution is indicated on 
the plot. 

mass planets could even be preferentially found orbiting metal- 
poorer stars. 

Using the numbers presented in Table|4] we measure a pre- 
liminary value for the Jupiter-to-Neptune ratio as a function of 
stellar metallicity. We estimate the value using the metallicity 
distribution of stars with only jovian planets in our HARPS 
sample, and compared it with the same distribution for stars or- 
bited only by Neptune-like planets. For this, we used the sam- 
ple of three stars mentioned above (HD4308, HD69830 and 
HD 160691, belonging to the sample presented in this paper - 
case 1), all stars only with neptunian planets listed in Table|4| 
(case 2), and all stars in the table (case 3). 

The results show that in the same three metallicity inter- 
vals mentioned above, the Jupiter-to-Neptune ratios are, by or- 
der of increasing metallicity, 5:1, 28:1, and 30:1 (case 1), 5:3, 
28:3, 30:0 (case 2), and 5:3 28:4 30:4 (case 3). In other words, 
this result tentatively suggests that the Jupiter-to-Neptune ra- 
tio may be higher at higher metallicities. This is expected by 
models o f planet formation based on the core-accretion process 
(iBenz et al. 2006). 

These final results should be interpreted with caution, how- 
ever, since most of the stars that have only with Neptune-like 
planets listed in Table|4]are M-dwarfs. It has been proposed that 
M-dwarf stars may have a lower rate of Jupiter-like planets, but 
(at least in the short planet period range) a high incide nce of 
Neptu ne-like planet detections (Endl et al. 2006; Bonfils et al. 
2007h . The above mentioned bias in the case of a search for 
neptunian planets around metal-rich (or planet-host) stars may 
also induce important biases in this analysis. 

8. Summary 

We have measured accurate stellar parameters for a well- 
defined sample of solar-type stars, which were originally ob- 
served as a part of a search for very low-mass planets. We have 
presented a catalog of stellar parameters for 45 1 stars which 



is available online. We compared our results with other works 
in the literature. The effective temperatures obtained were also 
compared with values derived using IRFM. All comparisons in- 
dicated a very good agreement in despite of the different meth- 
ods used by the different authors. 

A new calibration was developed for the determination of 
effective temperature using the index color B -V and metallic- 
ity. 

The results presented here were explored to study the 
metallicity-planet correlation, in particular for stars hosting 
planets in the Neptune-mass range. We showed that, when 
compared with their higher-mass counterparts, neptunian plan- 
ets did not follow the same trend. Low-mass planets ap- 
peared to be formed in a lower metal content environment. 
Interestingly, the Jupiter-to-Neptune ratio was found to be an 
increasing function of stellar metallicity. The results presented 
in this paper have shown that the study of the chemical abun- 
dances of Neptune-host stars provides new important con- 
straints on models of planet formation and evolution. 
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